Transport properties in the d-density wave state: Wiedemann-Franz law 
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We study the Wiedemann-Franz (WF) law in the d-density wave (DDW) model. Even though 
the opening of the DDW gap (Wo) profoundly modifies the electronic density of states and makes 
it dependent on energy, the value of the WF ratio at zero temperature (T = 0) remains unchanged. 
However, neither electrical nor thermal conductivity display universal behavior. For finite temper- 
ature, with T greater than the value of the impurity scattering rate at zero frequency 7(0) i.e. 
7(0) < T <C Wo, the usual WF ratio is obtained only in the weak scattering limit. For strong 
scattering there are large violations of the WF law. 

PACS numbers: 74.25.Fy, 74.20.Fg, 74.20.De 



In a recent paper Hill et a/.0 have observed large viola- 
tion of the Wiedemann-Franz (WF) law in (Pr,Ce) 2 Cu02 
driven into the normal state through the application of 
a 13 Tesla magnetic field. At very low temperature 
T < 0.2K, the thermal conductivity is found to be much 
less than the value estimated from the D.C. conductiv- 
ity. Above 0.3K the opposite holds. This observation 
suggests that an exotic state of matter may exist in the 
normal state of (Pr,Ce)2Cu02- Hill et al. consider spin- 
charge separation as one possibility. 

Recently (i-density jKasitijiDDW) order has received 
considerable attentioru'lJatffl as a possible exotic state 
of matter with a pseudogap which breaks time rever- 
sal symmetry because it introduces bond current with 
attendant small orbital magnetic moments. The pseu- 
dogap has d-wave symmetry. This is the symmetry ob- 
served in studies of the variation of the leading edge of 
the electron spectral ^density by angle-resolved photoe- 
mission spectroscopy,!!] as a function of angle in the Bril- 
louin zone in the normal state of underdoped cuprates. 
A pseudogap with d-wave symmetry implies important 
energy dependence of the quasiparticles density of states 
(DOS) at the Fermi surface (FS). Energy dependence in 
the DOS leads to impurity scattering rates that also de- 
pend on energy and the applicability of the usual WF 
law is no longer guaranteed. 

In this paper we consider the WF law within the DDW 
model. As yet, this model has not been shown to apply to 
the pseudogap regime of the cuprates. Here we take the 
point of view that nevertheless it can serve to understand, 
in this concrete case, how energy dependence in the DOS 
can alter the WF law. 

In the DDW state, the gap with magnitude Wq has 
d-wave symmetry and opens up at the antifcrromagnetic 
Brillouin zone of the Cu02 plane. Away from half filling, 
in the underdoped regime, the FS falls at the chemical 
potential fx (which would be zero at half filling) and we 
assume that \ii\ <C Wq. Provided that the effective im- 
purity scattering rate and temperature arer-also small as 
compared with Wo, a nodal approximations can be used 
to describe the electric as well as the thermal conductiv- 
ity. 

We consider a tight binding energy dispersion as 



a function of momentum k of the form: ek = 
— 2io[cos(fc a: ) + cos(fcj,)], where to is the in-plane hop- 
ping amplitude. At half filling the FS coincides with 
the antifcrromagnetic boundary where the DDW gap 
W"k = (W /2) {cos(k x ) —cos(ky)] opens up with amplitude 
W). Most properties of the DDW state are determined by 
the nesting vector Q = (it, 7T|i_.for example ek+Q = — £k 
and Wk+Q = — Wk- See Ref.Ertl for detailed properties. 
We begin with the Hamiltonian 



H = 



J dx^( x ) (_|_^ a ( x ) 

+ \J dxdy^(x)^(y)y(x-y)^(y)Va(x),(l) 

where ^( x ) creates an electron of spin a at x and 
y(x — y) is the electron-electron interaction. A spin sum- 
mation is implied. Using the definition of the DDW gap 
in momentum space: 

iW k = -^ k - k '(C k , +Q >0Ck'a("')} (2) 
k' 

one can obtain the mean field Hamiltonian of the DDW 
state. 

Let us consider the real part of the electrical conduc- 
tivity <r(Q). In the long wavelength limit (q — > 0) the 
current operator in momentum space j e (0, 51) is 



j e (0,fi) = -e^VkC+HfaCk^ + fi) , 



(3) 



where v k = de k /dk and C£(o;) = ^C kT (w), C k+QT (w) 
Note that we use four vector notation: k = (k, to) and 
J2k = EkEu- From the current-current correlation 
n(0, ifi), we have <r(fi) = — ^ImH re t(f2), where 

n(0, ift) = e 2 ^ vlTi \&(k, iCo + iQ)f 3 G(k, iw)f 3 ] (4) 

k 

Here the matrix Green's function G(k, ioj) is 

G(k,iuj = -2 (5) 

[iuj + n) 2 - E£ 



2 



with id) — iuj — Sq(«cD) and £k = \] e\ + IF, 2 . It 
has been assumed that the other components 
the self energy can be absorbed into ek and Wk 
It is useful to introduce the spectral functions 
Aij(k,u>) = — 2ImGy(k, ui+i5), for example, in(k,u) — 

i^+ll^i^r] + («k->«k,£flk-»--Ek), and 

Ai 2 (k,w) = 
where 7(w) 



f 2 7 (w) 
" W k«k | [( w + M +Sk) ! »+7(«) ! 'j 

= -ImE 0l ret(iw), U k = 



Note the differences between G sc for the DSC and G 
for the DDW. Using the spectral function A(k, ui) = 
— ImG(k, oj + iS) and B(k, to) = — InxF(k, u> + iS), where 
F(k, u>) is the anomalous Green's function, the D.C. con- 
ductivity becomes 

( flc^-iSk)} , *s C (T,0)=e*J2vl [A(k,.) 2 +i?(k,.) 2 ] 



|(l + e k /£ k ), 

and Vk = — e^/Ek). Now we obtain the D.C. 

conductivity cr(T, O = 0) as 

*(T,0) = e^l /^(-|£) [4u(k )W ) 2 -| J 4 12 (k jW )| 2 ] 
k ^ 

(6) 

where f(u>) is the Fermi function. In the nodal ap- 
proximation £ k ^ ^j/fif. £ k = pcos(6) and 
W k = psin(0). Then at T = we obtain a (0,0) = 
fi (%) A(0) , where 



A(u>) 



7(w) 



w + fl 



arctan 



7(w) 



(7) 



tr(0,0) depends only on 7(0) because only the zero fre- 
quency limit of A(w) enters at T — 0. 

This result shows that a(0, 0) depends not only on the 
chemical potential fi (and so on the filling) but also on 
the scattering rate 7(0). This is to be contrasted with 
the well-known universal value of the DC conductivity 
for the DSC: er sc (0, 0) = 2^(vf /v SCt9 ), where t> SCjS is the 
DSC gap velocity. For the DDW case a universal value is 
obtained only in the case when \i — > 0, which corresponds 
to half filling. In this limit cr(0, 0) reduces precisely to 
a sc (0,0) for the DSC with Wq playing the role of DSC gap 
(Ao). We see that it is because the DDW gap develops 
at the antiferromagnctic boundary rather than at the FS 
which is shifted by the chemical potential, which leads to 
the absence of universal behavior. 

It is instructive to contrast the DDW case with the 
DSC case in a more formal way. The charge current has 
the form for the DSC 

jSo(o, fi) = -e ]T VkV>kM4(^ + n) (8) 



where fp^ui) = (C kT (cj), C_ k |M). This leads to the 
current-current correlation 

n sc (0, ifl) = e 2 ^ vlTr [G sc (k, iu> + iO)G sc (k, iui) 

k 

(9) 

which is to be contrast with Eq. (|J). The matrix Green's 
function is: 

A /, ^tq + A k ri + (ek - /x)f 3 

G sc (k,uj) = - -2- (10) 

(iuj 2 - e k - A*) - A k 



(11) 



At T = 0, we obtain a sc (0, 0) = £ (^) A sc (0) , where 



A sc (u) = 2 



1 



arctan ■ 



7H 



j(u) 



(12) 



Next we consider the case of finite T in the range 
7(0) < T -C Wo- In this case id) ~ iuj — So(iw), 
namely, w can be vised in the evaluation of So to a 
good approximations Then £ ,ret(^) = J^q-i , where 

Go = ^ Ek with N ° bein S the D0 ° S at the 

FS, r is a scattering rate proportional to the impurity 
concentration, and c is the inverse of the impurity po- 
tential. For the Born limit c > 1 while in the unitary 
limit c — > 0. Applying the nodal approximation, one ob- 
tains Go — 



TT 2 N l 



-*f(a,+ M ) + ( W + M )ln(M) 



Thus for the Born limit we get 7(0;) = 70 t^ 1 , where 



7o 
7u 



r/c . For the unitary limit we have instead j(u>) = 
where « ir 2 T/4. These re- 



( Wo \ 



hr 



Wp 



suits for 7(0;) parallel the well-known results for the DSC, 
which are recovered when [i = with the DDW gap re- 
placed with the DSC gap. The most important feature 
of impurity scattering for our consideration of transport 
properties is that j(u>) acquires a frequency dependence 
and this leads to a violation with T of the WF law as we 
will see soon. For the Born limit a remarkable simplifica- 
tion for A(oS) occurs; namely, A(ui) becomes independent 
of frequency and this leads directly to no violation of the 

WF law. With A(lo) = 1 + + J^j arc t a n (^) ~ 



f because 70 <C Wq and it follows immediately that 



o-(0,T< Wo) 



e 

2^ 



Wo 

7o 



(13) 



which is temperature- independent. For the DSC in the 
same limit a sc (0,T < A ) ~ (e 2 /n)(v f /v sc , g )(A /j ). 
It is larger than the DDW results by a factor of two 
if Wq — Ao. The difference is traced to the fact that 
A sc (v) — 2 [(7r/2)A /7o]. This serves to illustrate that 
DSC and DDW order do not generally give the same an- 
swers. This is expected since in one case there is Cooper 
pair condensation while in the other there is none. 

We next consider heat transport in the DDW state 
since the WF law is a statement about the ratio of the 
thermal to electrical conductivity. The heat current j' i (x) 
can be calculated from the continuity equation: 7i(x) + 
V • j' l (x) = 0, where is 7i is the Hamiltonian density of 



3 



Eq. (Q). Define j h = } } j + jg, one can show that 



Thus the heat current-current correlation is 



and 



+ 53(^+2") uf C)9 Tr [G 8C (k,iw + ifi)fi& c (k ) icD)fi] (22) 

When the spectral functions are introduced, we arrive 
at 



V-#(x) = 9 / rfy^(y-x) ft«{x)1>Uv)Mv)Mx) 



^t(x)ipl(y)^/3(y)tp a {x) + h.c. 



(15) 



In momentum space, as q — > 0, 



1 ^ fcLu 2 ( 8f\ 

T2l^ J ^ {-Q^J 



j)(0, fl) = J] f w + - J v k C k (a;)faC7 k (a; + fi) (16) 

Applying mean field theory and keeping terms rel- 
evant only to the DDW order, we find iq ■ 
j£(q,n) = |(* a + A* g -- F 9 - Y* q ) where 

= Y,k ujW k+gC : l +cla (uj)C k+cla (uj + Q) and = 
A g (VKfc+q — > Wq) with a definition of the DDW gap 
Eq. (0). Now we obtain 



Applying the nodal approximation to Eq.(|20|), we ob- 
tain 



{vj + vlj [A(k,c) 2 -5(k,c) 2 ](23) 



7r 2 \_ 9w y l v fl u/ 



•AM (24) 



and exactly the same result holds for the DSC with ji — > 



and u Q 



As T -> 0, «/T = 



.4(0) for 



jg(0, n) = -*53 ( w + ? VffC£+q(w)f5,(7k(w + 0) the DDW while k sc /T = | f-^- + ^) for the DSC 



(17) 

where v g = dWk/dk. Therefore, the heat current be- 
comes 



j h (0,O)= J2 



n 



vkC k (^)T 3 c k ( w + n) 



2 
3 

In this case n sc /T is universal and does not depend on 
impurity scattering. In contrast n/T for the DDW has 
a dependence on 7(0) as well as on doping through the 
chemical potential. However, for the Lorenz number L — 
k(T)/[Tct(0, 0)] the scattering rate drops out and we find 



Ln = 



-iVgCt + Q(w)T 3 C' k (w + fl) 



3 c- 



1 + 



This shows that the WF law is 



( 18 ) obeyed at T = in the DDW state and the DSC case 



Note that we assume dv g /dt = so that we neglect the 
extra terms which depend on the time derivative of the 
gap velocity. 

The thermal conductivity /c(q = 0, 17) follows from the 
Kubo formula for the heat current-current correlation: 
= ~ WU^ m ^ret(^) > where 



R K (0,iQ.) =5Z( W+ 2) «kTr [<$(k,iw + in)f3G(k,w)f 3 ] 

k 

+ ^2{ UJ+ 2J [G(k, + in)f 3 G(k + Q, itD)f 3 ] (19) 



Again making use of the spectral functions, we obtain 



"(0) 
T 



-y 



jn 2 



du> 

L 

2tt 



duj 



[A u (k,u,) 2 -|A 12 (k, W )| 2 ] 



(20) 



Before proceeding further it is of interest to contrast our 
DDW derivations with the DSC case. For the DSC the 
heat current is 



jt(0,Q)= ^ 



" + 2 



(v g —> v sc ^ g ) and its value differs from the conventional 
one only by a very small correction of order (v g /vf) 2 due 
to a d-wave symmetry of the gap. 

A very similar result can be obtained in the case 7(0) < 
T <C Wo- In this regime we have already seen A(uj) ~ 
zm> for thc Born limit so OP) = 1 ( 2i + 2l) ffii. For 

the DSC, Wo — * A and w g — > «sc, ff - Because of Eq.(|l3|) 
for the DDW the Lorenz number reduces to the conven- 
tional value: L = ir 2 /(3e 2 ) for 7(0) < T < W . But for 
the DSC we arrive instead at the remarkable result that 
L sc = 7r 2 /(6e 2 ), a reduction of a factor of two. While 
we obtain this results analytically, Graf et aZJl3 have cal- 
culated L sc numerically and their work serves as a nu- 
merical verification of our result. Since the temperature 
scale for which this happens is 7(0) <C T, in the clean 
limit this switch-over from Lq to Lq/2 can happen at ex- 
tremely low T. In sharp contrast with the DSC, in the 
DDW case there is no change in the Lorenz number in 
the Born limit. 

It is not possible to obtain analytic results for the uni- 
tary limit. In general thc Lorenz number L(T)/L(0) is 
written as 



V k ^ k (w)T 3 'i/>k(w + 0) 



-Vg, SC ^ k (w)fl?/> k (w + ft) 



(21) 



L(T)_ 3 /Mf) 2 (-§^MH 
La tt 2 $du (_^ A{uj) 



(25) 
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FIG. 1: The normalized Lorenz number as a function of 
temperature (T). The dash-dotted curve is the result for 
the Born limit. Other curves are for the unitary limit. For 
the upper solid curve -yu /Wo = 0.001 and n/Wo = —0.01. 
The dashed curve is for the same chemical potential but for 
7c/ /Wo = 0.01. The lower solid curve is for 7(7 /Wo = 0.001 
but now n/Wo = —0.15. 



As we mentioned earlier, in the unitary limit "/(lu) = 
7c/ m ~ 2 ( ) ' Numerical results are presented 

in Fig. 1. We show results for four different cases which 
serve to illustrate what is possible. The upper solid curve 
is for Ju/Wq — 0.001 and /J. /Wo = -0.01 which shows 
a large peak around T/T* — 0.05. We have taken T* 
to be given by its mean field value: Wq/T* — 2.14 as 
in the DSC case. A very large positive violation of the 
WF law is seen. We need to point out, however, that 
while we have not shown a(T) and n(T)/T individually, 
in this case they both show large variations with T re- 
flecting the important frequency variation of 7(0;) for the 
unitary limit, which is not compensated for by the ex- 
plicit variation of A(uj) in Eq.(^). For the Born limit 
an exact compensation takes place so that A(lu) turns 
out to be a constant. This leads to the usual WF law 
with no T dependence, which is shown as a dash-dotted 



line in Fig. 1. For the dashed curve 7(7 /Wo = 0.01 and 
li /Wo = —0.01. Increasing 7(7 makes the deviations from 
the conventional Lorenz number smaller. The same ef- 
fect is obtained when \fj,\ is increased, effectively pushing 
the FS further away from the zero in DOS. The second 
solid curve has 7j//Wo = 0.001 but now h/Wq = —0.15, 
away from half filling. Now the deviation from the con- 
ventional Lorenz number can be negative as well as pos- 
itive depending on T but the amplitude of the violation 
is small because the DDW gap becomes less effective at 
changing the DOS near the FS. (Note that < W for 
the validity of the nodal approximation.) 

Our main conclusions are as follows. At T = 0, only 
the zero frequency limit of the imaginary part of the self- 
energy enters into the calculation of the electrical and 
thermal conductivity and the conventional Wicdcmann- 
Franz (WF) law is recovered. In contrast with what is 
found for a d-wave superconductor (DSC), for a d-density 
wave (DDW) state, neither electrical nor thermal conduc- 
tivity show universal behavior. Each depends on the im- 
purity scattering rate. But this dependence is the same 
and cancels from the Lorenz number as T — > 0. We 
were also able to obtain analytic results for low but finite 
temperature. In this case we found no change in the WF 
law for the Born limit even though the Lorenz number 
is reduced by a factor of two from its conventional value 
for the DSC. For the unitary limit, however, the Lorenz 
number increases rapidly at low temperature on a scale 
set by the zero scattering rate 7(0). In a case considered 
it rises above 2.5 around T/T* ~ 0.05 and then acquires 
a more moderate temperature variation. This case cor- 
responds to the chemical potential (/z) small compared 
to the DDW gap. When \/j,\ is increased sufficiently, the 
Lorenz number becomes approximately equal to its con- 
ventional value and its temperature dependence is small. 
It is important to realized that when the Lorenz num- 
ber is found to vary significantly with temperature, so 
do both electrical and thermal conductivities. This is 
generic to the model in which quasiparticlcs arc respon- 
sible for the. transport. Such a model cannot explain 
experimental in which the D.C. conductivity is almost 
independent of temperature while the Lorenz number is 
strongly dependent on it. 
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